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Hanbury Brown and Twiss (HBT) correlations, 
i.e. correlations in far-field intensity fluctuations, 
yield fundamental information on the quantum 
statistics of light sources, as highlighted after the 
discovery of photon bunching [l -3j. Drawing on 
the analogy between photons and atoms, similar 
measurements have been performed for matter- 
wave sources, probing density fluctuations of ex- 
panding ultracold Bose gases [4-8j. Here we use 
two-point density correlations to study how co- 
herence is gradually established when crossing 
the Bose-Einstein condensation (BEC) threshold. 
Our experiments reveal a persistent multimode 
character of the emerging matter-wave as seen 
in the non-trivial spatial shape of the correlation 
functions for all probed source geometries from 
nearly isotropic to quasi-one-dimensional (quasi- 
ID), and for all probed temperatures. The quali- 
tative features of our observations are captured 
by ideal Bose gas theory [9j, the quantitative 
differences illustrate the role of particle interac- 
tions. 

HBT correlations can be related to a quantum in- 
terference effect reflecting the multimode nature of the 
source [2 . In analogy to thermal light sources, atom 
bunching in expanding thermal Bose gases has been ob- 
served [IHH] . Well above condensation threshold thermal 
Bose gases are sufficiently dilute so that atom-atom in- 
teractions are negligible and ideal gas theory provides an 
accurate description. 

Extending the atom-photon analogy into the quantum 
degenerate regime, the absence of bunching has been ob- 
served in an out-coupled atom-laser jTOl [TT] and in an 
expanding BEC ^ . This suggested a perfect coher- 
ence of these systems, similar to a monomode optical 
laser where interferences leading to HBT correlations are 
essentially absent, yielding only spatially and temporally 
uncorrelated poissonian shot noise [12]. 

Yet, studies of first order coherence properties of Bose 
gases near the BEC phase transition have revelled the 
importance of thermal excitations, reducing the coher- 
ence also below threshold [131 El- Moreover, atom- 
atom interactions have been identified as a cause for 
a multimode nature in very elongated degenerate Bose 
gases. Originally predicted for the limit of weakly in- 
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FIG. 1. Fluorescence imaging of density correlations, a, 
Schematic of the experimental geometry, b, Example of the 
density distribution (lateral cut, log scale) of a Bose gas 
slightly below the Bose-Einstein condensation threshold af- 
ter 46 ms expansion time, c, Profile of b, integrated along y. 
The red line indicates a fit to the density distribution of the 
thermal fraction, yielding the temperature T. 



teracting degenerate ID systems [15], similar behaviour 
is also present for very elongated 3D degenerate Bose 
gases [16l [17]. This effect has been demonstrated ex- 
perimentally through measurements of density fluctua- 
tions [18] or two-point correlations [19] for temperatures 
significantly below the BEC threshold. 

In the following, we probe the second order correla- 
tion function g2 of expanding Bose gases across the BEC 
phase transition, at zero and finite distances, to map 
out the gradual establishment of matter-wave coherence. 
Varying the source geometry, we explore the regime from 
'3D' to 'quasi-lD' physics. Close to the BEC threshold, 
the fluctuations of the most populated quantum modes 
are significant and common theoretical models describing 
interacting Bose gases, as in [15-.i-17j. fail. We therefore 
restrict our analysis to a qualitative comparison to ideal 
Bose gas theory predictions, allowing us to cover all in- 
vestigated regimes. 
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FIG. 2. Density correlations results. a,b,c, Density correlation functions of expanding Bose gases n2 (lateral cut in {Ax, Ay) 
plane) at 0.99 T^, 0.95 T^, and 0.85 T^, the aspect ratio A of the atomic source is 13.5. cl,e,f, (circles) Radial cuts of a,b and 
c. (dashed line) Radial cuts of the corresponding autocorrelation of the mean density profile, (solid line) Radial cuts of ideal 
Bose gas theory's predictions for the second order correlation function ^2(0, A^). g,h,f, (circles) Radially averaged axial cuts 
of a,b and c over 160 /xm where the shot noise peak 5 is excluded. The width of the exclusion region is 32 /xm. (dashed line) 
Radially averaged axial cuts of the corresponding autocorrelation of the mean density profile, (solid line) Radially averaged 
axial cuts of ideal Bose gas theory's predictions for the second order correlation function Q2{Ax). 



The function ^2(1*, r') measures the probability of joint 
detection of two particles at positions r and and hence 
relates to the density fluctuations of the system and their 
spatial correlations. For an ideal Bose gas above the BEG 
threshold ^2(0) = 2, highlighting an excess of density 
fluctuations (bunching). For finite distances, ^2(1*, r') de- 
cays to unity on the length scale of the temperature- 
dependent coherence length. For a true monomode 
source ^2(1*, r') = 1, demonstrating perfect coherence. 
In general the shape of ^2(1*, r') reflects the interplay be- 
tween the occupied modes of the system through the spa- 
tial scales of density fluctuations. 

To obtain the second order correlation function of a 
Bose gas after its release from the trap, we use our novel 
fluorescence imaging [20|. Its high spatial resolution, 
single-atom sensitivity and exquisite signal-to-noise ra- 
tio enables us to probe the ^2 function with an accuracy 
at the percent level. We record a thin slice of the atomic 
density in the horizontal x-y plane at the central part 
of the gas after 46 ms expansion (see Fig. [T^-b). The 
thickness of the slice is set to 225 /im by adjusting the 
duration of the excitation pulse to 500 /is. 

Our matter- wave source is a ^'^Rb gas of adjustable 
temperature prepared in the |F, m^) = 1) state in 
an atom chip magnetic trap [21j. The chip design allows 
us to control the parameters of the magnetic trap and 



hence the shape of the source over a large range of aspect 
ratios A = ojy^z/^x from 4 to 120, where ujx and ujy^z are 
the axial and radial trap frequencies respectively [22^. 

Modifying the aspect ratio of the trap we explore dif- 
ferent regimes of Bose gases (see supplementary informa- 
tion). For the most isotropic traps (A = 4.75, 13.5), gases 
significantly below the BEG threshold can be considered 
in the Thomas-Fermi regime with a chemical potential ja 
of the order of 3 to 4 huy^z- In this case the expansion 
starts with a hydrodynamic phase. For the most elon- 
gated trap (A = 113), gases with temperatures below the 
degeneracy threshold have ja slightly smaller than hwy^z, 
entering the weakly interacting quasi- ID regime [23j. In 
this case the expansion of the system is essentially colli- 
sionless [17 . 

We obtain the temperature of the system T and de- 
duce the number of atoms in the thermal fraction of the 
gas A/th by fitting the thermal tails of the axial {x axis) 
density profiles using an analytical formula based on ideal 
Bose gas theory (see Fig.jl]::). By comparing the fit result 
to the total observed signal, this also yields the number of 
condensed atoms (if present) within the measured slice. 
A simple model for the expansion of the degenerate frac- 
tion of the gas yields A^bec, the total number of condensed 
atoms. Finally, for each set of temperature T and total 
atom number N = A/th + A/'bec, we estimate the criti- 
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FIG. 3. Density correlations comparison. a,b,c, Peak height of the second order correlation function of expanding Bose gases 
for various ratios T/T^ and three different trap aspect ratios A = 4.75 (light blue), 13.5 (orange) and 113 (green). It is defined 
as 1x2(0) (see caption Fig. [2]). The solid lines indicate the predictions obtained from ideal Bose gas theory. The light gray area 
display the range where the source is fully thermal. Note that the total atom number N is not held constant between the 
data points, however this is taken into account in the value of . d,e,f, (circles) Radially averaged axial cuts of the density 
correlation function n2 (see caption Fig. [2]) which parameters correspond to the colored square displayed respectively in a,b 
and c. (dashed line) Radially averaged axial cuts of the corresponding autocorrelation of the mean density profile, (solid line) 
Radially averaged axial cuts of ideal Bose gas theory's predictions for the second order correlation function 02 (Ax). 



cal temperature of the system through the relation 
= a^^^Tcj where Tc{T^N) corresponds to ideal Bose 
gas theory predictions for the critical temperature. The 
experimentally obtained factor a depends exclusively on 
the source aspect ratio A (see Methods). 

Averaging over typically one hundred experimental 
repetitions with identical experimental parameters, we 
calculate the mean autocorrelation of the density profiles 
and normalize it by the autocorrelation of the mean den- 
sity profile. We obtain the normalized density correla- 
tion function 1x2, which can be decomposed into the sum 
of two terms, 5 and 02- The term ^2 contains the desired 
information about two-particle correlations whereas 5 is 
a contribution due to atomic shot noise (see Methods). 

Typical results of density correlation functions of ex- 
panding Bose gases from a moderately anisotropic source 
(A = 13.5) are shown in Fig. |2] (the corresponding graphs 
for a nearly isotropic trap (A = 4.75, Fig. SI) and a 
quasi-lD gas (A = 113, Fig. S2) can be found in the 
supplementary information). The atomic shot noise 5 
appears as a dark oval spot at the center of the function. 
Its eigen-axes are rotated by 45° with respect to the axes 
of the atomic source x, corresponding to the direction 
of the imaging laser beams. The structure corresponding 
to 02 is visible behind this central spot. Its anisotropy 
reflects the aspect ratio of the source gas. 

For a thermal gas above the BEC threshold (T > T^) 
we observe bunching. In this region, where inter-particle 
interactions are weak, we find an almost perfect agree- 



ment between the experimental observations and the sec- 
ond order correlation function ^2 simulated within ideal 
Bose gas theory. As input parameters we take the ther- 
modynamical properties obtained from the data used to 
compute n2 and account for the imaging resolution |9|. 

Below the critical temperature (T < T^) the main ex- 
perimental observations valid for all explored aspect ra- 
tios (see Fig. [2] and Figs SI and S2 in the supplementary 
information) are: 

(1) The establishment of coherence along the radial di- 
rection: The RMS width of the Q2 peak along the y axis 
rapidly grows when crossing the critical temperature 
until it saturates due the finite radial size of the system 
(see Fig. [2]i-f). This can be seen as an experimental ob- 
servation of the establishment of radial coherence [14^, 
where the shape of n2 changes from a decreasing expo- 
nential (Fig. |2]i-e) to a profile set by the spatial shape of 
the condensed cloud after expansion (Fig.[2]f). The RMS 
width of the Q2 peak is related to the radial coherence 
length in the trapped system through the hydrodynamic 
expansion of the gas, which imposes a scaling factor on 
the radial size of the cloud (see Methods). 

(2) 02(0) > 1 and the appearance of a dip below unity at 
finite distances along the axial 02- Due to the anisotropy 
of the trap, thermal excitations of the system will typ- 
ically be spread over more modes axially {x axis) than 
radially. Below but close to the BEC threshold, few of 
the lowest lying axial modes will be macroscopically oc- 
cupied and the shape of the second order correlation func- 
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tion 02 (Ax > 0) can then be interpreted as a result of the 
interference of all contributing modes. Hence the time-of- 
flight expansion implements a very sensitive heterodyne 
detection of weakly occupied modes of the matter-wave 
source. The dip below unity at finite distance Ax along 
the Q2 axial profiles is a direct consequence of this inter- 
ference (see Fig.|2ji-i). This observation is similar to ob- 
servations reported in recent experiments and theoretical 
work on weakly interacting quasi- ID Bose gases [171 HI] 
(for a comparison see Fig. [sjd-f). Most interestingly we 
find this behaviour - generally associated with quasi- ID 
physics - also for the most isotropic trap probed (aspect 
ratio A = 4.75). 

We would like to point out that (2) corrects the widely 
established image of a perfectly 'fiat' correlation func- 
tion for a Bose gas immediately below the BEC thresh- 
old El [ini [n]. Such behaviour highlights the infiu- 
ence of the non-ground-state modes (thermal depletion) 
of the system, whose population saturates below the BEC 
threshold. Working with small atom numbers and our 
highly sensitive detector allow us to accurately probe the 
transition regime and the graduate saturation of excita- 
tions when crossing the BEC threshold. 

To give an overview we display the behaviour of n2 (0) , 
set by 02(0), for all three measured aspect ratios and 
temperatures in a wide range of different ratios T/T^ in 
Fig. |3^-c. Most strikingly we find 02(0) > 1 for all tem- 
peratures and aspect ratios, which indicates a persistent 
multimode nature even for a three-dimensional BEC far 
below the condensation threshold. Lowering the tem- 
perature we observe a slow decrease of the amplitude of 
HBT correlations, which illustrates the gradual reduction 
of the thermal depletion of the system. 

Ideal Bose gas theory reproduces the basic physics and 



the qualitative features of the experimental observation 
for moderate aspect ratios as discussed above also be- 
low the critical temperature. We attribute the remaining 
quantitative deviations between this description and our 
observations to particle interactions (see also Fig. |3]l-f 
and Figs SI and S2 in the supplementary information). 

For highly anisotropic sources with a large aspect ratio 
(A = 113) close to a quasi- ID system {ja < hw)^ we ob- 
serve no reduction of 52(0) with temperature (see Fig.jsJ: 
and Fig. S2 in the supplementary information). This be- 
haviour is expected for very elongated systems where low 
lying axial excitations of the system remain macroscop- 
ically occupied also significantly below the degeneracy 
temperature p!5UT7l fT9] . It is based on the dominant 
infiuence of particle interactions, hence ideal Bose gas 
theory fails to describe even the qualitative behaviour at 
T < 

Our experimental findings call for a more complete 
theoretical description of interacting Bose gases at the 
threshold to Bose-Einstein condensation, where fiuctu- 
ations of competing modes are significant. With such a 
theory at hand, measurements of density correlations will 
allow a detailed quantitative characterization of the mode 
occupation of the source. For equilibrium systems well 
below the BEC threshold this can be used for precision 
thermometry, where methods based on the observation 
of a thermal background fail (as demonstrated for deeply 
degenerate ID systems [HI US])- For ultralow tempera- 
tures well below the chemical potential, this would allow 
to probe quantum fiuctuations, e.g. quantum depletion, 
of 3D or ID Bose gases. Our studies can directly be ex- 
tended to non-equilibrium systems, where fundamental 
questions on equilibration, thermalization and integra- 
bility arise |24H26]. 
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creases as well. Here Tc — huk^^ (A^/C(3))^^^ where N is 
the total number of atoms, uJ = (cJccCJ^cj^z)"*^^^, h and ks are 
Planck and Bolzmann constants respectively and (^{x) is the 
Riemann zeta function. To ensure thermal equilibrium of the 
gas we keep a constant RF knife for 200 ms at the end of the 
evaporative cooling ramp. 



Comparison with ideal Bose gas theory 

The constant a introduced in the main text is adjusted such 
that cooling the gas to = a^^^Tc coincides with the experi- 
mental observation of the appearance of a degenerate fraction 
of atoms (A^bec ^ 0). The value of a is slightly smaller than 
unity for the moderate trap aspect ratios (A = 4.75,13.5), 
which qualitatively agrees with theoretical studies of the shift 
of the critical temperature of interacting Bose gases [27| . For 
the most elongated trap (A = 113), a is of the order of 3, 
which qualitatively agrees with predictions for the degener- 
acy temperature for ID systems [28] . 



Fit of the axial profile of the thermal fraction 

Using an analytic expression for the density flux / of a 
freely expanding ideal Bose gas (see [9 , last formulas on page 
3 and 5) and assuming an elliptical Gaussian shape for the 
excitation beams, hence neglecting any saturation effect, we 
deduce a formula describing the axial profile of density " slices" 
of the thermal fraction of expanding Bose gases 
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27TW 



dydz exp 
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that we can use to fit our experimental results. Here to is the 
expansion time, zq the position of the centre of the excitation 
laser beam profiles, 2 At the exposure time and w the RMS 
width of the excitation beams. For the data presented in this 
work to 46 ms, 2 At 500 /us and = 10 /xm [20]. This 
model depends on two parameters, the temperature of the gas 
T and its fugacity, which fixes the total number of atoms in 
the thermal fraction of the gas Nth- To avoid the divergence 
of the population of the ground state of the system when the 
fugacity reaches its maximum, only the excited states of the 
ideal Bose gas are considered for the calculation of pis{x). 



METHODS 

Experimental Setup 

We prepare ultracold gases of a few 10^ ^^Rb atoms in the 
|F, ttif) = |1,— 1) state confined in one out of 3 different loffe- 
Pritchard type magnetic trap obtained on top of a multilayer 
atom chip The 3 harmonic trap parameters are uox = 

27V X 20 Hz axially and cj^,^ = 27r x 2260 Hz radially (A =113), 
Ua: = 27T X 160 Hz and ujy,^ = 27r x 2180 Hz (A = 13.5) and 
uja: = 27T X 320 Hz and ujy,z = 27r x 1520 Hz (A = 4.75) 
respectively. 

We use forced radio-frequency (RF) evaporation to cool the 
gas close to or below the Bose-Einstein condensation thresh- 
old. As a direct consequence, colder gases contain less atoms 
and their corresponding critical temperature = a'^^^Tc de- 



Estimation of A^bec 

For traps with moderate aspect ratios (A = 4.75, 13.5), the 
condition /a ^ huOx,y,z is fulfilled quickly after passing the 
BEG threshold (see Fig. S3 in the supplementary information) 
and the condensed part of the system can then be assumed 
to be in the Thomas Fermi regime. Scaling laws can then 
be used to describe the gas expansion ^29l [30] . We find scal- 
ing factors of 6 (axial) and 625 (radial) for A = 13.5 and 28 
(axial) and 430 (radial) for A = 4.75 respectively. With such 
model, deducing the number of condensed atoms A^bec from 
the number of condensed atoms counted in a density slice is 
straightforward. 

For the most anisotropic trap (A = 113), we use an isotropic 
Gaussian ansatz for the radial density distribution of the de- 
generate fraction of the gas. Fitting the RMS width of the 
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Gaussian along the y axis, it is then straightforward to infer 

iVbec 



Density correlation function 

The density correlation function is the sum of the two con- 
tributions 02 and 5. The first term 02 can be expressed as 
02 = ^ * o^, where ^ (Ar) is an average of the second order 
correlation function ^'2(1*, r^) over the spatially inhomogeneous 
density profile of the gas with Ar = r — r^ kept fixed, is 
the effective two-particle point spread function (PSF) of the 
fluorescence imaging and * is the convolution product. The 
function accounts for the diffusion of single atoms during 
the imaging process and optics aberrations, which blur the 
position of single atoms in the images 20 . The other con- 
tribution to the density correlation function 5 is due to the 
atomic shot noise and is proportional to the inverse of the 
mean density of the gas. 

Here we slightly extend the analytic expression of the sec- 
ond order correlation function of an expanding ideal Bose gas 
derived in 9 (see top right formula on page 6) to take into 
account the imaging scheme. Using a 2D isotropic Gaussian 
model for we obtain an analytic expression for 02- 



SUPPLEMENTARY INFORMATION 

Second order correlation results for A = 4.75 and 
A = 113 

Here we present figures showing the same analysis data 
as for Fig. [2] in the main text but for different trap aspect 
ratios: A = 4.75 for Fig. [4] and A = 113 for Fig. jt] 
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FIG. 4. Density correlations results. a,b,c, Density correlation function of expanding Bose gases n2 (lateral cut in (Ax, Ay) 
plane) at 1.05 T^, 0.94 T^, and 0.84 T^, the aspect ratio of the atomic source A is 4.75. d,e,f, (circles) Radial cuts of a,b and c. 
(dashed line) Radial cuts of the corresponding autocorrelation of the mean density profile, (solid line) Radial cuts of ideal Bose 
gas theory's predictions for the second order correlation function 02(0, Ay). g,h,f, (circles) Radially averaged axial cuts of a,b 
and c over 160 /xm where the shot noise peak 5 is spatially excluded. The width of the exclusion region is 32 /xm. (dashed line) 
Radially averaged axial cuts of the corresponding autocorrelation of the mean density profile, (solid line) Radially averaged 
axial cuts of ideal Bose gas theory's predictions for the second order correlation function 02 (Ax). 
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FIG. 5. Density correlations results. a,b,c, Density correlation function of expanding Bose gases n2 (lateral cut in {Ax, Ay) 
plane) at 1.00 T^, 0.95 T^, and 0.85 T^, the aspect ratio of the atomic source A is 113. cl,e,f, (circles) Radial cuts of a,b and c. 
(dashed line) Radial cuts of the corresponding autocorrelation of the mean density profile, (solid line) Radial cuts of ideal Bose 
gas theory's predictions for the second order correlation function 02(0, A^). g,h,f, (circles) Radially averaged axial cuts of a,b 
and c over 160 /xm where the shot noise peak 5 is spatially excluded. The width of the exclusion region is 32 /xm. (dashed line) 
Radially averaged axial cuts of the corresponding autocorrelation of the mean density profile, (solid line) Radially averaged 
axial cuts of ideal Bose gas theory's predictions for the second order correlation function 02 (Ax). 



Thermodynamical properties of the data 

In Fig. [6] we show the thermodynamical properties of 
the different data used to compute the second order cor- 
relation functions. As explained in the methods section, 
the total atom number N and the number of condensed 
atoms A^bec are obtained from the individual images using 
simple models for the expansion of the gas. According 
to ideal Bose gas theory, the condensation threshold is 
characterized at a given temperature by a critical atom 
number 




over which the excited states of the system are saturated 
and the atoms accumulate in the ground state of the 
trap. This prediction has to be rescaled using the phe- 
nomenological factor a to comply with our experimen- 
tal observations, accounting for the effects of interactions 
and dimensionality. 

We estimate our error on the temperature estimation 
of the mean density to be less than 10 %. This gives an 
error on the estimation of the ratio T/T^ of the order 
of 10 %. Nevertheless, close to T^, the measurement is 
much more accurate since the presence of a few condensed 
atoms is easily visible in individual images. This feature 
is used for the determination of a. 

To estimate /i in the case of the most isotropic trap 
(A = 4.75, 13.5) we assume the system to be in the 
Thomas Fermi regime, in which case 



/iTF 



15 



1 2/5 



where a = 5.3 nm is the s-wave scattering length and 
uj = {uJx'^yUJzy^^ . Such assumption is only strictly valid 
when /iTF ^ ^i, with i = x^y^ z, which is fulfilled in our 
case only for the coldest data. 

In the case of the most elongated trap, we assume the 
system to be in the quasi- ID regime and estimate ja fol- 
lowing Ref. [23 . Here again, the approach is only strictly 
valid for the coldest samples. 



10 




-I 1 1 1- -I 1 O— I- -I 1 OEOCCtO- 

0.5 1 1.5 1 2 0.5 1 

T [/iK] T [/iK] T [fiK] 



FIG. 6. Thermo dynamical properties of the data. a,b,c Scatter plot of the means of the derived total atom number N (grey 
disks) and total condensed atom number A^bec (red disks) versus the temperature T, following the procedure described in the 
Methods section. The error bars give the standard deviation of N and A^bec over individual samples. The dashed line indicate 
the ideal Bose gas theory prediction for the critical atom number (see text). The solid line gives the same prediction corrected 
by the phenomenological factor a. d,e,f Reproduction of a,b,c zooming on the behaviour of A^bec- g,h,i Scatter plot of the 
estimate of the chemical potential of the gases. The color of each disk is related to the ratio Nhec/N. The method used to 
obtain fi is described in the text. It is strictly valid only for the darkest disks. The dashed line indicates the value of the radial 
trapping frequency Uy^zj"^^- 
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Normalisation of the fluctuations 

When averaging the fluorescence pictures in order to 

obtain the density correlation function n2, we average 

clouds with slightly different total atom numbers. These 

fluctuations can be corrected by normalizing 1x2 by the 
2 

factor 1 + Var(S')/S' , where S is the total fluorescence 
signal in each picture, S its mean value for the considered 
clouds and Var(S') its variance. This factor differs from 
one by less than a percent for all the data presented here. 
No further normalization is necessary in order to obtain 
n2 equal to unity at large distances Ax, Ay. 



Imaging resolution 

While the resolution of the fluorescence imaging along 
the z axis is fixed by the waist of the light sheet beams 
and the duration of the pulse, in the x — y plane it is to a 
large extent determined by the diffusion of atoms during 
the detection process Hence, each individual atom 
appears as a fluorescence pattern^ composed of typically 
5 detected photons randomly distributed over a few ad- 



jacent pixels. The broadening of the atomic shot noise 
peak 5 is given by the mean autocorrelation of the fluo- 
rescence patterns. Its elliptic shape reflects the increased 
diffusion of the atoms along the direction of the excitation 
beams. It is modelled by a 2D Gaussian with eigen-axes 
rotated by 45° and fitted RMS widths V2 x 5.5 /im and 
/im. However, 5 is insensitive to the random off- 
set of each fluorescence pattern centroid with respect to 
the actual position of the corresponding atom. This off- 
set has to be included in a full description of the imaging 
resolution, yielding the two-particle PSF o^. 

Assuming a 2D isotropic Gaussian model for o^, we can 
estimate its RMS width from the observed axial width of 
02 above the condensation threshold. The RMS width of 
related to the axial coherence length of the gas [6 , 
is expected to be much smaller than the observed RMS 
width of 02- Thus, we assume the latter to originate 
only from the finite resolution (see Fig. [2^). This, to- 
gether with the finite resolution of the fluorescence imag- 
ing along the z axis, explains the reduction of 52(0) below 
2 for systems close or above the BEG threshold [6 . As 
the observed width stems from a two-particle correlation, 
it is sensitive to the actual resolution of the system, in- 
cluding the centroid offset. Finally, we obtain a value of 
V2 X 8 /im for the RMS width of 



